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Systems of differential equations with state-dependent delay are considered. The delay dynamically 
. . depends on the state i.e. is governed by an additional differential equation. By applying the time 

04 , transformations we arrive to constant delay systems and compare the asymptotic properties of the original 

' and transformed systems. 
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Delay differential equations successfully used to model and study a number of applied problems in 
physics, biology, chemistry. Basic mathematical theory of constant delay equations can be found in the 
classical monographs [51 [5] and references therein. 

Differential equations with state-dependent delay (SDD) attracted much attention during last decades 
and many deep results have been obtained so far (see |9l [11] [151 [10] and references therein). Such equations 
with discrete state-dependent delays are always nonlinear by their nature. As described in [9], this type 
of delay brings additional difficulties in proving such basic properties of solutions as uniqueness and 
' continuous dependence on initial data. The main approach to get the well-posed initial- value problem is 

to restrict the set of initial functions and hence the set of solutions to -functions [5]. For an alternative 
approach where an additional condition on the SDD is used to get a well-posed initial- value problem in the 
space of continuous functions see |121 113| . In this note we rely on the classical approach [5] and compare 
the SDD problem with another one with a constant delay. This constant delay problem is constructed by 
using the so-called time transformations [21 [3] . This transformation could be applied to any particular 
solution along which the deviating argument is monotone. To be assured that the monotonicity holds for 
^ I all solutions we concentrate on the system when the state-dependent delay is governed by an additional 

CO ' differential equation and provide a sufficient condition for the monotonicity of the deviating argument. 

I This type of equations is used to describe some models of population dynamics, see [1] and references 

therein. In [T one could also find motivations to study this type of SDD and comparison with the 
frequently used case when SDD is presented by explicit or implicit functionals. 
■ Our main goal in this note is to compare the asymptotic properties of the SDD system with the 

corresponding ones of the system after the time transformation. We try to find conditions which guarantee 
that such properties as stability, boundedness and compactness of the initial SDD problem survives under 
the time transformations. 

2. Time transformations 
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We study the following non- autonomous system with state-dependent delay (see the autonomous case 



yit) = f{t,y{t),y{t-r^m, t > t\ (1) 
m = -y^{ri{t)-^ + G{y{t)), t>t\ (2) 



with initial data 



y{t) = g{t), te[t"-h,t% (3) 
Vit') - ry". (4) 

Here y G i?™, rj G R, rf , rj > 0. The function 77 is a state-dependent delay since it is a solution of the 
equation ([2]) where there is a dependence on y. 

In the sequel we will denote h = 2rj and also X = C^{[—h, 0]; i?™) x R with the natural norm. 

Lemma 1. Let f and G he Lipschitz and \G{y)\ < firj for all y £ i?'". Then for any g G C^{[t^ — 
h, t'^]; R™), r](t'^) = rf G [0, 2rj\ the system (Op-Q^ has a unique global solution such that ri{t) G [0, 2?)] for 
all t > t'^ . The solution continuously depends on initial data {g]rf). 
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Proof of lemma 1. The existence of solutions is simple since the righthand sides of equations ((T]), ^ 
are continuous. Solutions are global due to Lipschitz properties of / and G. The uniqueness follows from 
the well-known results on the state-dependent delay equations (see e.g. [9]) since we consider Lipschitz 
initial function g. 

Using the property |G(y)| < /i^, one can easily show that (for any y) any solution of ([2]), Q satisfies 
T](t) G [0,2^ provided T](t°) = E [0,2?j\. 

Now we show the continuous dependence on initial data. For the simplicity of presentation we put 
t° = 0. Let us consider a pair {g; 77°) e C^{[—h, 0]; i?™) x [0, h] C X and an arbitrary sequence (5"; ?7°'") 
such that ||(g";r/"'") - {g;Tl°)\\x -> as n ^ 00. 

We rewrite the system (UJ-Q in the integral form 

y^t) = g"(0) + f f{s, y^s), y-{s - fj" {s)) ds, (5) 
Jo 

77"(i) e-^*(77°'" -?7) + / e-^(*-'')G'(2/"(s))ds. (6) 

"'0 

Similar equations are for initial data {g; ff^). For the differences of solutions, using the Lipschitz properties 
of / and G (the corresponding Lipschitz constants are denoted by and Lq) we have 

l2/"(i)-y(i)|< 1.9" (0) -3(0)1 

+Lf f {\y-{s - ry«(s)) - y{s - 77"(s))| + \y{s - r;"(s)) - y{s - f]{s))\ + \y-{s) - y{s)\} ds, 
Jo 

|77"(i) - m < e-^'\v"'- ~ f)°\ + Lg f |y"(s) - y{s)\ ds. 

Jo 

Let us fix any T > 0. Since all solutions are in time (see e.g. [9])) for i > we can denote by Lg^T the 
Lipschitz constant of the solution y{t),t £ [—h,T]. Hence \y{s — ri^{s)) — y{s — fj{s))\ < Ly,T\'rf'{s) — fj{s)\ 
for all s e [0,t] C [0,T]. Denoting for short ^"(i) = max^g[o,t]{|y"(T) - {/(r)] + |77"(r) - f]{T)\} and 
C"^ = max{2L/ -I- Lq] Ly^r} we obtain 



0</3"(i)</3"(0)+i/ inax |g"(r)-g(r)KC^ / /3"(s)ds. 
Te|-/i,oj 

We apply the Gronwall's inequality and get for all t £ [0,T] 

max { |y" (r) -y(T)K by" (t)-^(t) I } 

TG|0,tJ 



exp(t • max{2L/ -f Lg; Ly^r})- 



< \il + Lf) max |g"(r)-5(r)| + hO'"-7t 

The last convergence and equations ([l}, ([2]) give a similar estimate (convergence) for the time derivatives, 
so we finally get 

I2/" - y|ci[o,t] + I*?" - jylcMo^t] -> as n^+00. 
This gives the continuous dependence on initial data and completes the proof of lemma 1. ■ 

For any solution (y; 77) of the system ([Ij-lll]) we call the function a given by 

a{t) = t-ri{t), t>t" (7) 

the deviating argument of {y',ri). 

Our goal is to investigate properties connected to the time transformation approach introduced in 
[U [3]. We are going to use a function t = a{s) called time transformation [3] to convert a particular 
solution {y; 77) of the system (HJ-Q into a solution (z; x', ct) of constant delay system 

i(s) = f{a{s),z{s),z{s — h))a{s), s > s°, 

z{s) = g{io{s)), se[s"~h,s% 

xis) - -^l{xis)^v)ais) + Gizis))a{s), 

^ X(s°) = 77°, 
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where a satisfies the algebraic equation 

/ ol{s) - x(s) = a(s - /i), s > s°, , , 

Here cj : [s° — ft,, s°] — s> i? is an arbitrary C^-function with positive derivative and such that a;(s° — h) = 

Remark. We notice that equation 0) is different from the corresponding rules used in and 
since here we have no a given lag function. 

The time transformation a is constructed step by step (see ([T])) by the rule 

a{s) = a'\a{s - h)), s e [s° + (k ~ l)h, s° + kh], /c = 0,1,2.... (10) 

Here we used (see ([7])) a{a{s)) — a{s) — r]{a{s)) — a{s) — x(s) = a{s — h), since x{s) — ri{a{s)). 

It is clear that one needs the invertibility of a to define a. In [3] the rule was used assuming that 
&{t) > (or <j{t) < 0). More precisely, it was used for those solutions along which &{t) > (or &{t) < 0). 
In our study we can give a simple condition which guarantees that along all solutions we have &{t) > 
and hence a is always invertible. Such a simple condition is < 1. It is easy to see using ^ that in 
this case |77(t)| < n\ri — rj\ + \G{y{t))\ < firj + /iry < 1. Here we used the assumption |G'(y)| < ^rj. Now ([7]) 
implies &{t) = 1 — ri{t) > 0. 

Remark. Taking into account that the state- dependent delay ?/ takes values in [0,2rj\ — [0,h], one 
could say that the assumption < 1 means the "slow changing of delay" in the range [0, 2rf\. 

It is important that if a{t) > 0, then d(s) > since on the initial time segment a{s) — uj{s) > 0, s £ 
[fiO - h, s"] and a{s) = ^^i"'^^' = a{a{l-h)) ^(^ ^ ^) ^ > step by step (see also [3, p.28]). 

Here we used (j{a{s)) = a{s — h) (see ([5]), (dU])). 

By construction (see [21 propositions 1,2]) the connection between a solution {y,rj) of ©-(H) and the 
corresponding solution (z,x, a) of ([S])-® is given by 

y{t) - z(a-i(t)), t e [tO-/i, +oo), z{s) = y{a{s)), x{s) = vHs)), s e [s°-ft, +oo), t = a{s), = a(s°). 

(11) 

Lemma 2. Let f and G be as in lemma 1. Consider a sequence {(<?"; f?*^'")} such that \ \{g^\rf'^) — 
(5; Tf)\\c'^{[t'>-h,t%R"^)>^R ~^ sequence {w"} such that | (I>| |,7i[(0_/j ^o] as n — !■ 00. Then for 

any S > Q the sequence of time transformations a" uniformly converges to a i.e. max^g^sO 30^5] |a"(s) — 
a(s)| — > as n — > cxD 

Proof of lemma 2. First, using we have from lemma 1 that Taajii^^t" ,t° +t\ |o'"(^) ~ as 

n 00. This convergence and ^ imply that ||cr" — CT||(7i[to_to_|_T] 0. 

Since -^^(t) > on [i", t° + T], then there is a ^ > such that ^a{t) >2S>0 for aU t e [t°, t'^ + T]. 
This and Her" — o'\\ci[tQ,to+T] — ^ give 

— CT"(t) > (5 > for alH e + T] and n > n^. (12) 

dt 

Using the definition of a (see (|10p ) and the convergence w" ^> w as n — >■ 00 we only need to show the 
property max.re[o,T] I(o'")~^(t) — as n — >• 00. Let us denote by 7"(s) = (ct"')~^(s),7(s) = 

(a)^^(s). Assume opposite i.e. 7"(s) 'p, 7(5) on some [s°,s'' + S\. Hence 

3eo>0:ViVeN 3nAr>7V 3 s„„ e [5°, 5° + 5] : |7"(s„„) - 7(s„„)| > Eq. 

Considering = 1, 2, ... we get two sequences {nfej^i and {s„^ C [s°, s° + S\ such that 

l7"Hs„J-7(s„J| >£o. (13) 

Since \s^\ s^ + S\ is compact we have s € [s°, s'^ + S\ and a subsequence again denoted by {f^fe}fe^Ll such 
that s„j. — !• s e [s", s*^ + S"]. We can write 

t"^^™.) - 7(s«J = (7"'=(snJ - 7"'=(s)) + (7"'= (5) - 7(s)) + (7(s) - 7(snJ) 
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The last term vanishes due the continuity of 7, the second one due to the point- wise convergence 7"'(s) — >■ 
7(5) for all s G [s°, + S]. Hence the only possibility to satisfy is that there is an integer ki such 
that for all k > ki one has 

\r'M-rHs)\>eo/2. 
The last property together with s and differentiability of all 7"'' imply that the derivatives ^7"*° 

are unbounded in a neighborhood of s. This contradicts the property ^cr"(t) > 5 > (see ((H])) since 
^7"(s) = TT^T^- The proof of lemma 2 is complete. ■ 

Now combining lemmata 1 and 2 (and property ijrj < i) we can formulate the first result on the 
continuous dependence of the time transformation on initial data. 

Theorem 1. Let f and G he Lipschitz and \G{y)\ < ^rj < ^ for all y £ i?'". Consider a sequence 
{((7"; ?7°^")} such that 1|(5"; — {g',v'^)\\c^{[t'>-h,t°];R'")xR ~^ a sequence {w"} such that ||a;" — 
to\\c^[t'>-h,t"] — > as n — 00. Then the time transformation gives the sequence of the corresponding 
solutions {(z",x",a")} of the constant delay system (^-(^ (see Ull) ) such that for any S >0 one has 

max {\\z-{s) - -z{s)\\ + |x"(s) - x{s)\ + |a"(5) - a{s)\} ^ as n ^ 00. 

se[sO,sO+S] 

and a"(s) > for s e [s", s° + S]. 

Remark. We should notice that (0) is the system of coupled differential and algebraic equa- 
tions. It is necessary to comment on how to solve it. The way is different from the one 0/ since 
we have no a given lag function (c.f. |3 Section 2.1]). Using Q), we write for s € [0,/i] (and than 
continue step by step with the step h): Qf(s) = x(s) + w(s — h). Then we substitute it into the differential 
equation for x m |3i to get x{s) = —IJ'ixi^) ~ v) (x(*) + ^{^ ^ h)) + G{z{s)) (x{s) + i1j{s — h)). Hence 
x{s) [1 + /i(x(s) G'(2:(s))] = {^m(x(s) ~ ^t) ) + G(z(s)} cj(s — h). We remind that the assumption 

2^rj < 1 implies |/i(x(s) —rj) — G{z{s))\ < 1. It gives 

m - {-A^(x(s) G{z{s)} co{s ~h)[l + ^l{x{s) -V)- G{z{s))]-' . (14) 

Now to rewrite the first equation in (0) we use again a{s) — x{s) +'^{s — h) (and a{s) — x{s) +Lb{s — h) ), 
substitute it in Q) and use |J^[ ). It gives z{s) — f{a{s), z{s), z{s — h)) a{s) — /(x(s)+Cli(s — /i), z{s), z{s — 
h)) ixis) + Lois - h)) ^ fixis) +c.{s- h),z{s), z{s - h)) + 1) . _ /,) ^ /(^(,) + 

Lo{s - h), z{s), z(s - h)) [1 + m(x(s) - ^ - G{z{s))]-' ■ uj{s - h). 
Hence can be rewritten for s € [0,h] (the first step) as 

' z{s) = fixis) + u{s - h), z{s), zis ^ h)) [1 + ^lixis) ~ v) - G{z{s))]-^ -Lois-h), se[0,h], 

z(s) = .g('^(s)), se[0,h], 

X(s) = {-f4xis)-7r))+Gizis)}ujis-h)[l + f,ixis)-7i)-Gizis))]-\ 

[ x(o) = v", 

(15) 

It easy to see that the last system gives solution (z(s),x(s)) for s S [0, /i]. Then the time transformation 
a is found by a{s) = xi-^) + ^{s — h). In general it reads as 

iis) = fixis) + ais-h),z{s),zis-h))[l + tiixis)-v)-Gizis))]-^ -ais-h), s > s° , 
z{s) = giujis)), se [sO-/i,sO], 

X(s) = {-fiixis) -v)) + G(z(s)} [1 + m(x(s) -v)- Gizis))]-' a(s - h), s> s°, 
X(s°) - 7f, 

(16) 

and can be solved step by step. 

3. Connection between asymptotic properties of systems (Il|)-(l4]) and ([8|)- (j9l) . 

In this section we discuss how to determine if some qualitative properties of solutions of the initial 
state-dependent delay system ([l])-(|4l) survive the time transformation i.e. still be valid for the correspond- 
ing solutions of constant delay system ([S])-®. We are also interested to connect the known properties of 
solutions of dH])-® with the ones of ©-dH). 
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Let us start with the discussion of the property of the exponential stabihty. For the simpHcity of 
presentation we assume that the zero function y[t) = is a solution of Hence, by PT|) . z(s) = 

will be also a solution of ([HI)-®- 

Assume one has 

Ms)\\ < i?"e-^'(^-^")||z,o||c([-„,o];fl"),s > s",D^,D' > 0. 
Hence, by (jlip . we get 

It is easy to see that if (and only if) " is bounded, then we have 

llyWII <i^'e-^'(*-*")||z,o||c([-M];fl'")- 

The above considerations show that the exponential stability of the zero solution z{s) =0 implies the 
exponential stability of the zero solution y{t) = provided there are positive constants D^,D^ such that 
D'^t - D^a~^{t) < for all t > t°. Similar estimates give the inverse implication i.e. the exponential 
stability of the zero solution y(t) = implies the exponential stability of the zero solution z{s) = 
provided there are positive constants C^,C^ such that C^s — C^a{s) < for all s > s°. Since by ((TT|) 
t = a(s), we arrive to the following 

Definition 1. We say that "s-time" and "t-time" are equivalent if there are constants > 0,A^ > 
0, e R,B'^ e R such that A^t + < s < A^t + B^. 

Remark A. It is evident that in this case we also have -^s — ^ l£ t < -^s — ^ . 

We saw that in the case of the equivalent "s-time" and "t-time" the property of the exponential 
stability survives under the time transformation. Another consequence of the equivalence is that t — ^ -|-oo 
if and only if s — > +oo, which is clearly important for the study of long-time asymptotic behavior of 
solutions. 

The last result suggests to study the notion of time-equivalence in detail. Let us try to find if in our 
case we have the equivalence. The rule ([TOl) and ^ show that we need to analyse the function a^^. Using 
the property a{t) > t—h (bounded delay) and invertibility of a, we get ^^^{t) < t+H. Hence, by PU)) . one 
has a{s) = a~^{a{s-h)) < a{s-h) + h. Particularly, a{h) < a{0) + h,a{2h) < a{h) + h < a{0) + 2h,etc. 
Hence the property a{kh) < a(0) -I- kh and the strict increase of a give the following estimate 

a{s) < {a{0) + h) + s. (17) 

Since a{s) = t, the estimate (jl7p means the lower bound in definition 1 and the upper bound in remark A 
(with A^ = 1, Bi = -(a(0) -I- h)). 

The supplemented bounds in definition 1 and remark A are less obvious. For the moment we do not 
claim that it is true in general case, but present an additional assumption which guarantees the bounds. 

Let us assume that the value of delay is bounded from below by a positive constant, say hi > 0. 
More precisely, rj{t) > ft-i G (0,^ for all t > Q. A sufficient condition for the last property is |G(y)| < 
IJ.\'q — hi\ for all y G i?™. Under the above condition we have (T^^(t) > t + hi. By (|10p . one has 
a(s) = a^^{a{s — h)) > a{s — h) + hi. 

Particularly, a{h) > a(0) -I- hi,a{2h) < a{h) + hi < a{0) + 2hi,etc. Hence the property a{kh) < 
a(0) + khi and the strict increase of a give the following estimate 

a{s)>{a{0)-hi) + ^s. (18) 

Combining ([T7)) and ([T5)) we conclude that "s-time" and "t-time" are equivalent (with A^ = 1,B^ = 
-(a(0) + h),A'^ ^ Ji-^B^ ^ -■^(a(O) - hi) in definition 1 and remark A). 

Having the equivalence proved we can use it to compare the asymptotic behavior of the corresponding 
dynamical systems (process), constructed by solutions of systems before and after the time transforma- 
tions. 
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Let us consider an autonomous case of ([T]) i.e. the system (c.f. 



(19) 
(20) 



with initial data 



yit) 

77(0) 



g{t), te[-h,0], 



(21) 
(22) 



We can restrict our study (using lemma 1) to the set of initial data 



Xf^ = {(<?,77°)|<?(0) = r(g(0),g(-r;0))} cCi([-^0];i?™) ^ [o^;^]_ 



(23) 



The set Xfa is an analogous to the solution manifold used in [TS] (see also [9J). We notice that the 
reason for this restriction is to have smoothness of solution at zero i.e. y(0— ) = 5(0—) — j)(0+) = 
/°((;(0), (7(— 77")). It is easy to see that Xfo^ is invariant. 

We define the evolution operator S°'{t) : Xfa — > X/a, associated to the system P^ -(P ^ . by the 
formula S°'{t){g]rf^) = {yt;ri{t)), where (y; 77) is the unique solution of P^ - (P^ . It is easy to see that 
under our assumptions the pair {S°',Xfa.) constitutes a dynamical system (in other words, the IVP 
p^ - (P^ is well-posed in Xfa). For more definitions and details on dynamical systems see e.g. [SlfT^H], 

To discuss the properties of solutions to the non- autonomous system (|S]), let us remind the 
following definition from |16[ pp. 112-119]. Let i? be a Banach space. Consider the two-parameter family 
of maps {U{t,T)}, U{t,T) :£'—!> E, parameters t € R,t > t. 

Definition 2. [16, p. 113]. A family of maps {U(t,T)} is called a process on E if 

(i) U{t,t) — I =indentity; 

(a) U{t, s) o U{s, t) = U{t, t) for allt>s>T e R. 

Since ^ has constant delay h we have no need to restrict our study to a solution manifold. We 
define 



and define U{s,s'^){(p;r]'^;Lu) — (z^; x(s); o^s), where (2;;x;a) is the unique solution of ([8]), ([9]) with initial 
data (iy9; 77*^; oj), here ZgO = ip. 

We notice that to come back to the system ([T^ - (P^ . using a particular solution of ([HI, (HI), one 
restores function g as follows g{t) — Zso{a~^{t)) for t ^ [t'^ — h,t'^]. 

Let us continue to discuss which asymptotic properties of (IT])-(|3]) survive the time transformation i.e. 
still be valid for the corresponding solutions of constant delay system ([5])-([n]). 

One of the important properties of dynamical systems and processes are boundedness of solutions and 
their compactness (asymptotic compactness) see e.g. [8l[14l|4]. Below we always assume that t — >■ -|-cxd if 
and only if s — ^ +00, which is true, for example, in the case of the equivalence oit—time and s—time. One 
can easily see, by (jlip . that ||j/(i)|| < C,t > ti is equivalent to t|2(s)]| < C,s > si. Hence, the existence 
of a bounded absorbing set for Zg -coordinate is equivalent to the existence of a bounded absorbing set 
for yt -coordinate, both in the space C{[—h, 0]; i?™). To go further, let us discuss the following additional 
assumptions on the time transformation a: 

(Al) 3Ci'" > 0,Vs > si =^ d(s) < C^'". 

(Al') 3C2'° > 0,Vt > ii ^d-i(i) < C^'". 

(A2) a is uniformly continuous on [si, +00). 

(A2') a^^ is uniformly continuous on [ti, +00). 

(A3) a is uniformly continuous on [si, +00). 

(A3') is uniformly continuous on [ti, -l-oo). 

Using (fTTj) . we have z(s) = y(a(s))d(s). Hence, the existence of a bounded absorbing set for the 
yt -coordinate in the space C^{[—h,0];R"^) implies the existence of a bounded absorbing set for the Zg 
-coordinate, provided (Al) is satisfied. The inverse implication is valid provided {Al') is satisfied. For 



£; = C([-/i,0];i?") X [0,h]x 
x{w e C^[~h, 0] I uj{-) > 0, uj{-h) = uj{0) - 77°; uj{0) [1 + ^(77" G{ip{0))] = io{-h)} 
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the system the existence of a bounded absorbing set means it is dissipative (for more details on 

this property see e.g. [Hlllllll])- Let us now assume that for t > ti the yt -coordinate belongs to a (pre-) 
compact set in the space C{[—h,0];R™). By the Arzela-Ascoli theorem theorem, the family {yt}t>ti 
is uniformly bounded and equicontinuous. Using (|lip. we have |^(s^) — z(s'^)\ — \y{a{s^)) — j/(a(s^))|. 
This and the above discussion show that for s > si the Zs -coordinate belongs to a (pre-) compact 
set in the space C{[~h,0]; R™), provided (A2) is satisfied. The inverse implication is valid provided 
{A2') is satisfied. The similar considerations in C^([— ft., 0]; i?™) need the estimate |i;(s^) — ^(s^)| = 
|y(a(si))d(si)-y(a(s2))d(52)| < \y{a{s^)) -y{a{s^))\\a{s^)\ +\y{a{s^))\ |d(si) - d(s2)|. We see that if 
for t > ti the yt -coordinate belongs to a (pre-) compact set in the space C^{[—h, 0]; -R™), then for s > si 
the Zs -coordinate belongs to a (pre-) compact set in the space C^([— /i, 0]; i?™), provided (A1)-(A3) are 
satisfied. The inverse implication is valid provided (^1') - (^3') are satisfied. Particularly, we shown 
that assumptions (A1)-(A3) and {Al') - {A3') connect asymptotic properties of the dynamical system 
{S''{t),Xfa) and the process C/(s,t) : E ^ E. 

Remarks. 1. Discussing assumptions (Al)-(A3), one could think that the family {as}s>si may 
belong to a (pre-) compact set in the space C[~h, 0] or even C^[~h, 0], but it is never true since a{s) — t 
is time, which is naturally unbounded. 

2. One can see that (Al) implies {A2), but {A2) ^ (Al), similarly {Al') => {A2'), but {A2') i> {A\'). 

3. We notice that {Al) gives a{s) < C^'"s + and similarly {Al') a^'^{t) < C^-"t + P. Using 
definition 1, we see that {A\), {AV) imply the equivalence of "s-time" and "t-time". 

Acknowledgments. The author thanks S. Maset for useful comments on an earlier version of the 
manuscript. 
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